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Many  different  types  and  configurations  of  nutation 
dampers  are  known,  and  a  single  theory  capable  of  describing 
the  action  of  all  of  them  would  be  very  complicated  indeed. 
However,  many  of  them  can  be  described  by  considering  the 
effect  of  a  small  off-axis  weight,  mounted  so  it  can  rotate 
around  the  axis  of  symmetry  of  the  body  on  a  bearing  with 
friction.  The  mercury  damper  can  be  described  in  this  way 
when  it  is  so  designed  that  the  mercury  does  not  spread  out 
in  the  channel  or  break  up  into  beads . 

In  the  derivations  which  follow,  the  following  postu¬ 
lates  are  assumed. 

(1)  The  body  carrying  the  damper  has  an  axis  of  sym¬ 
metry,  and  is  called  a  top  in  this  paper. 

(2)  The  damper  consists  of  a  small  point  mass  located 
a  distance  r  from  the  axis  of  the  top  and  displaced  a  dis¬ 
tance  t  along  the  axis  from  the  c.g.  of  the  top. 

(3)  The  damper  r.an  rotate  around  the  axis  of  the  top 
on  a  bearing  which  has  small  viscous  friction. 

(4)  No  external  moments  are  applied  to  the  body  or 
damper  mass . 


Coordinate  System 

In  order  that  the  time-varying  moment  of  inertia  terms 
be  eliminated,  the  coordinate  system  chosen  rotates  with 
the  damper  about  the  axis  of  the  top.  It  is  illustrated 
below: 

t 
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Fig.  1.  Moving  Coordinate  System 

The  origin  of  the  coordinate  system  is  the  c.g.  of 
the  top,  which  has  mass  M.  The  angular  velocity  components! 
of  the  top  are  and  u^,  while  the  angular  velocity  k 

components  of  the  coordinate  system  and  damper  are  cc,, 

and  ol,.  The  coordinates  of  the  mass  are  t,  0,  and  r,  while 

the  coordinates  of  the  c.g.  of  the  system  are  l,  0,  r 


The  Laeraneian 

Since  there  are  no  external  moments  and  the  bearing 
friction  is  not  conservative,  the  Lagrangian  of  the  system 
la  8 imply  the  rotational  kinetic  energy. 
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12  12 

T  «  T  4-  "o  Mvm  +  mv  ,  where: 
o  2  M  L  nr 

T  -  total  rotational  kinetic  energy 

Tq  ■  Aa?  4-  Bo^  4  j  Bu|  -  kinetic  energy  of  top 

VM  "  velocity  of  origin  of  coordinate  system 

vm  *  velocity  of  damper  mass. 

Define 


O  ■  OL  4-  oj  j  4  co  k  «  angular  velocity  of  coordinate 
^  system. 


rm  *  +  rk  “  position  of  damper. 


rcg  -  ^  -  position  of  c.g.  of  system. 


VM  ’  ’0x  rcg  ’  M+S  (0x  rm> 


vm  -  "x  <4  ’  rcg>  '  ^x?m> 

T  ■  ^  a4 +  \  B<^ +  4) +  i  «(?&$  <”  * 


-  ^  Ai^  +|!(i^+^)+|*(l>«  rm)*  , 


l 
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p 


where  m 


mM 

rn-t-M 


’’reduced”  mass  of  damper. 


O  x  rm  -  irc^  +  j*  -  rO)  +  k(-iup 
(Ox  r*m)2  -  r2^  +  i2(t^  +  (|)  +  r2rf  -  2trc^a 

T  “  I  A“k  +  1  +  f  (B  +  mi2)  (i^  +  t^) 

+ 

Define 

—  -  2 
B  ■  B  +  m4 

r  "2 

C  «  mr 

D  ■  mrt 

T*iA4  +  ic^+i  ®(u^  +  «ib  +  i  c4  '  "H0 

i 

Euler  Angles 

The  moving  coordinate  system  can  be  obtained  from  an 
inertial  coordinate  system  :coyozQ  by  three  rotations: 

(1)  An  angle  $  about  xQ 

(2)  An  angle  0  about  the  new  y 

(3)  An  angle  <p  about  the  new  x 
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The  additional  spin  of  the  top  can  be  expressed  in  terms 
of  an  additional  angle  17  about  x.  The  angles  are  illustrated 
below: 


*0 


Fig.  2.  The  Euler  Angles 


The  angular  velocities  are  given  in  terms  of  the  Euler 
Angles  and  17  by: 


0  -  +  $  cos  0 


UL  m  TJ  +  fl  -  »?+<p  +  ♦  COS  0 
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sin  8  sin  ip  +  6  cos  <P 


*  ip  sin  0  cos  <P 


The  equations  of  motion  are  obtained  from  the  kinetic 
energy  by 


u  ST  ST  M 
dt  Sqt  “  “  Mqi 


where  the  q’s  are  the  Euler  angles  0,  6,  and  <P  and  the 
additional  angle  rj.  The  M^'s  are  evaluated  by  the  principle 

of  virtual  work. 

In  the  case  where  there  are  no  external  moments,  all 
the  M^s  vanish  except  M^,  which,  for  viscous  friction,  can 

be  written 


M 


V 


-k q,  where  k  is  the  friction  torque  coefficient. 


The  differentiations  required  to  evaluate  the  equations 
of  motion  are  somewhat  complicated,  requiring  evaluation  of 
all  the  partial  derivatives  of  Cl,  ,  u  . 


The  results  of  this  process  are: 

(1)  The  rj  equation 

^  '  0) 

(2)  The  (p  equation 

+  C(  0  -  -  D(  +  0  uy)  -  o 
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(3)  The  0  equation 


[AcScak  +  B(^  "  A  ttjz)  +  c(wy  +  0  <*fe) 

+  »crf  -  «■£>]  cos  <p  -  ^-Auj^u^  +  B(o^  +  O  u^,) 
-  D(h  -  a^a^)J  sin  <p  «  o 


(4)  The  ip  equation 


cos  0  [ASc  +  C(h  -  a^o)z)  -  D(o^  +  O  wy)J 

+  sin  0  j^-Au^u^  +  3(0^+  O^) 

-  D(0  -  us  co  )\  cos  <p  +  /au>  us 
y  *  >  l  x  z 

+  B(t^  -  0(^)  +  C(u^,  +  O  u^.) 

4-  D(0^  -  u£)|  sin  <pj  -  o 


Fortunately,  it  is  possible  to  simplify  these  equations'^ 
considerably  by  combining  them  so  as  to  eliminate  the  trigo¬ 
nometric  functions  of  <p  and  0. 

Combining  (1)  and  (2), 


(5)  Au^  -  -k(u^  "  0) 


(6)  C(0  -  u^)  *  °(<*^  +  0  u^)  *  k(u^  -  0) 
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Combining  (2)  and  , 


[“AtScuV  +  +  A  wy)  "  D(0-  wyw7)J  cos  (p 

+  [AScWz  +  B(^  “  0  ^7)  +  C(a^  +  0  u^) 

+  D(0^  -  u?)J  sin  <p  -  o 


Combining  this  with  (3), 


(7)  Au^o^  +  ^  Wj)  +  C(u^  +  Q  a^)  +  D(f?  -  u£)  -  o 


(8)  -Awxu^  +  B(wz  +  O  (^)  -  D(Q  -  c^u>z)  -  o 


Solving  these  for  the  highest  derivatives, 


(9)  Sc  “  X(Sc  ‘  n) 

n  ~  “y“i  -  §(“i  +  0  V  +  c(“x  ‘  0) 

<*  +  0  “y  ■  i<o-  + 1  <s,^  • 

(10)  CJ-  _^k  ■;(%  -  (1)  +  s,“  +  W.K 

BC-D  x  BC-D  x  y  y  * 


.AC 


D k 


(11)  ul  -  “0  uj  +  --^“T  'jcui  +  -  -  0) 

y  BC-D  x  y  BC-i)  x 

U2)  vr£d“v-^-*i»i-r-(i)!-  <i) 

J  oxT  *  n  .r  A  *  air*  ** 


B4C 


B-tC 


B+C 
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For  a  small  damper  and  small  nutations, 

"z  St’  0 
C,D  «  E,  A,  B  -  B 

Under  these  conditions,  the  above  equations  reduce  to: 

(13) 

(14)  b  -  -  Q  -  o^)  + 

(15)  ^  -  O  |  - 

(16) 

The  effect  of  these  approximations  has  been  shown  to 
be  negligible  for  reasonable  dampers  by  analogue  computer 
simulations . 


T  r  an  s  f  Q  tma  t  Lon  _  to.  Pis  p_l  ac  gmgnlfi 

Equations  (13)  to  (16)  are  sufficient  to  define  the 
angular  velocities  l  it  it  is  desirable  to  determine  the 
angular  displacements. 


Set  X  •  j.in  6  cos  ^ 
Y  •*  sin  0  sin 


These  coordinates  are  rectangular  distances  between  the 
inertial  xq  axis  and  a  point  unit  distance  above  the  c.g.  of 

e 
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the  top,  on  the  axis  of  the  top.  If,  initially,  if:  -  <p  *  0  - 
o,  the  damper  is  on  a  line  parallel  to  the  -X  axis. 

Differentiating  these  expressions, 

*  •  • 

X  «  0  cos  0  cos  if)  -  if}  sin  0  sin  if) 

Y  »  0  .  os  0  sin  if)  4-  if)  sin  0  cos  if) 

4  •  •  * 

X  cos  if)  +  Y  sin  if)  *  0  cos  0-0  for  small  0 

•  4  • 

-X  sin  if)  +  Y  cos  if)  *  if)  sin  0 

Referring  to  the  Euler  angle  transformation,  the 
following  relations  are  evident: 

0  -  ojy  cos  <p  -  u>z  sin  <p 

if)  sin  0  *  co  sin  +  icz  cos  <p 


«  ■  <*  • 

0  -  (p  +  if)  cos  0  —  <p  +  if) 

•  • 

Combining  the  expressions  for  0  and  if)  sin  0, 
•  * 

X  cos  if)  +  Y  sin  if)  ■  oj  cos  <p  -  sin  p 


-X  sin  if)  +  Y  cos  0 


Co  sin  <p  +  co  sin 

y  * 


Solving  for  X  and  Y 


X  «  ujy  co3  («p  +  0)  -  «z  sin  (<p  +  (/») 
« 

Y  -  sin  (<p  +  0)  +  cos  (<p  +  if)) 
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Define 


(p  +  j/)  •  R 

f  -  R 


X  -  co  cos  R  -  co  sin  R 

y 


% 


sin  R  +  co  cos  F 

7. 


Similarly,  solving  for  c^  and 


‘■y 


X  cos  R  +  Y  sin  R 


*  "X  sin  R  +  Y  cos  R 


•  •  • 


“y 


-  ;x  +  R\)  cos  R  +  (Y  -  RX)  sin  R 


M  •  • 


U)  *  -(X  +  RY)  sin  R  +  (Y  +  RX)  cos  R 


Substituting  in  equations  (13)  -  (16)  yields,  after 
come  simplification: 


'•  v  *  An  • 

R  -  -  |(R  -  c^)  +  (X  cos  R  +  Y  sin  R) 

••  a  *  D  *  2 

X  +  g  o^Y  -  g  R  cos  R 

•*  a  •  n  4  2 

Y  -  f  uj^X  -  -  f  R  sin  R 
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2 

Set  Z  -  X  +  iY,  where  i  -  *1 


The  simplified  equations  of  motion  become: 


/-.-7N  n  k  ,  ,  AD  -iR  ,  -t  iRN 

(17)  R  »  -  ^  (R  -  u^)  +  -2gQ  u>x(Ze  +  Ze  ) 


.A  '  D’2  iS. 

(i8)  ^  '  ij  V  "  '  BR  e 


Solutions  of  the  Equations  of  Motion 

Equations  (17)  and  (18)  are  very  nonlinear,  so  that  a 
general  solution  cannot  be  obtained.  It  is  possible  to 
find  particular  solutions  or  approximate  solutions  which  will 
be  useful  in  describing  the  effects  of  the  damper. 

Three  cases  appear  to  be  useful: 

Case  I.  The  damper  does  not  move  relate  '  to  the  *~oy. 

The  solution  is  a  limit  cycle. 

Case  II.  The  damper  oscillates  with  small  amplitudes 
about  the  Case  I  solution.  This  is  called  the  ’’slow"  damping 
mode  in  this  report. 

Case  III.  The  damper  rotates  at  the  nutation  frequency. 
This  is  called  the  "fast”  damping  mode. 

Other  solutions  obtained  on  the  analogue  computer  appears 
to  be  mostly  mixtures  of  Case  'T  and  Case  III  types,  but 
one  other  is  of  considerable  interest,  This  is  extremely 
rapid  damping  which  can  occur  if  the  initial  excitation  is 
of  the  proper  size  and  phase  relative  to  the  damper.  This 
is  called  ''resonant”  damping  and  has  not  been  demonstrated 
analytically. 


Case  I:  Steady  State  Solution 
Choose  R  »  u^t 


Z 


-  ae 


1<4ct 


12 


V 


TN  4065-63 


Substitution  into  equation  (17)  and  (18)  shows  that  (17) 
is  satisfied  if  a  is  real  and  (18)  is  satisfied  if 


-D 
“  A-B 


the  solution  is 


R  ■  Sc* 


This  form  should  be  expected,  since  the  damper  behaves 
like  an  ordinary  unbalance. 


Case  II:  Oscillating  Solution 
Choose  R  -  +  a,  where  a  is  small 


Z  «  Z  +  Z.,  where  Z  is  the  solution  of  the 
o  1'  o 


homogeneous  equation  and  Z^  is  a  small  perturbation  on  this 
solution. 


z. 

« 

z 

1 

o 

1 

From  equation  (18) 


* 

Z 

o 


0 


2L  “  a  e 
o  o 


A 

t 


For  convenience, 


a0  is  chosen  to  be  real. 
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From  equation  (17) 


•  -5  “  ■  ^4a0  sin  ^bW 


-  4*^1 a  - 1 


D  20lSctn  .  .  v 
“  B  Sce  (1  +  ia) 


Solution  of  these  equations  is  straightforward  but 
somewhat  laborious.  The  results  are: 


a  »  a 


A2D 


o  BC(A-b)  Sc  [  ~1 


Ac  „„„  /A-  Bv  . 
C  cos  (“H4 


+  (^Sc  •*»  (T2)^] 


, ,A-Bv  _ 
.  »2n  „  1  B  ^'HcC 

7  *D  I.  iA  DaoSce 

zi  Tb  11  “ 


\  4C(A-B)  (2B-A)  -  i(MH) 


ScC 


.,,2  2  .A  . 

4cao _ ‘bV 

,k  .  ,  ,A-B,  te 


2(A-B)BC(£  +  K^jj2)^) 


r. 
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The  first  term  in  will  be  recognized  as  the  steady 

state  value  with  a  small  oscillation  superimposed.  This  term 
will  be  neglected,  since  it  does  not  influence  the  damping. 

Z  -  Z  +  Z, 
o  1 

..  L _ ^ _ l1!  V 

0  2(A-B)BC(£  + 


Analog  computer  solutions  of  equations  (17)  and  (18) 
indicate  that  tnis  estimate  of  tc  is  approximately  correct. 


15 


% 


TN  4065-63 


Case  III:  Nutation  Frequency  Solution 


A 

Choose  R  ■  -g  o^t  +  <p,  where  <p  is  a  constant. 

Z  *  Z  +  Z-.,  where  IzJ  «  jz  I 
01'  Ii  o 


a  e 
o 


.A  . 
i^oo  t 
B  x 


irom  equation  (17) 


-k,A-JL  ,  AD  .A  2  (  +i<p  ip, 

R  "  °  “  C  (  B  )  Sc  +  2BC  iB  We  ‘  e  ^ 


sin  „  .  MAiBIb 

A  DaoSc 


Since  sin  <p  *  1,  a  as  in  order  that  this  solution 


"°  A2D 


Sc 


exist. 

From  equation  (18) 


y  .Ai  D  /A  v2 /bS^  i<p  ,  .  DA2  2  i<p  4 

zi  ■  ti  "  ■  b  (b  e  e  *  p  s;e  e 


i-^co  t 
x 
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,A  *- 

Since  e  is  a  solution  of  the  homogeneous  equation, 

choose 


a^te 


.A  _ 
iB“xt 


Substituting  this  form  into  the  differential  equation; 


i-^  ^  (sin  cp  -  i  cos  <p) 


»  -DA 


+  i 


Z  *  Z  +  Z1 
o  i 


Z 


2k(A-B)c 

AB 


I 

L 


Since  aQ  is  not  a  common  factor  In  this  expression  no 

assumption  of  an  exponential  form  for  the  damping  is  reason¬ 
able.  Analog  computer  studies  of  the  equations  and  this 
solution  indicate  that  it  is  valid  for  a  short  time,  but 
that  a  much  better  approximation  is  obtained  by  "starting  * 
over"  after  each  nutation  cycle. 
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Define 


At^  «  time  per  nutation  cycle 


2v 


B 


ku> c 


ty  2  2 

I  t«o  o 


71 2  2 

Z I .  *  -  a„ 

I  t-Atn  o 


2k(A-B)Atn  ,DA  ,2  ? 

_ty  V  Atn 


AB 


zi2 


MUn  -Wl, 


At 


n 


2  ttDAu: 


-2k(A-B)  _ 

AB  +“7 


iz  2  a 


2  /2k(A-B) 

o  -  -  —p—r 


The  analog  computer  solutions  indicate  that  this  form 
is  reasonably  accurate  until  the  amplitude  has  decayed  to 
quite  near  the  critical  value. 
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